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We summarize recent results regarding the equilibrium and non-equilibrium 
behavior of cold dilute atomic gases in the limit in which the two body scatter- 
ing length a goes to infinity. In this limit the system is described by a Galilean 
invariant (non-relativistic) conformal field theory. We discuss the low energy 
effective lagrangian appropriate to the limit a — > oo, and compute low energy 
coeflicients using an e-expansion. We also show how to combine the effective 
lagrangian with kinetic theory in order to compute the shear viscosity, and 
compare the kinetic theory predictions to experimental results extracted from 
the damping of collective modes in trapped Fermi gases. 
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1. Introduction 

Over the last ten years there has been remarkable progress in the study of 
"designer fluids" , dilute, non-relativistic Bose and Fermi gases in which the 
scattering length between the Bosons or Fermions can be continuously ad- 
justed. In the following we are particularly interested in Fermi gases, since 
these systems arc stable for both positive and negative values of the scat- 
tering length, including the strongly correlated limit in which the scattering 
length is taken to infinity. 

The scattering length is controlled through a Feshbach resonance. Alkali 
atoms such as ^Li and ""^K have a single valence electron. When a dilute 
gas of atoms is cooled to very low temperatures, we can view the atoms as 
pointlike particles interacting via interatomic potentials which depend on 
the hypcrfinc quantum numbers. A Feshbach resonance arises if a molecular 
bound state in a "closed" hyperfine channel crosses near the threshold of an 
energetically lower "open" channel. Because the magnetic moments of the 
open and closed states are in general different, Feshbach resonances can be 
tuned using an applied magnetic field. At resonance the two-body scattering 
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length in the open channel diverges, and the cross section a is limited only 
by unitarity, (T{k) = ^Tr/k"^ for low momenta k. In the unitarity limit, details 
about the microscopic interaction are irrelevant, and the system displays 
universal properties. 

Near a Feshbach resonance the scattering length behaves as 



where ao is the non-resonant value of the scattering length (typically on the 
order of the effective range of the interatomic potential), B is the magnetic 
field, Bq the position of the resonance, and AB the width. A small negative 
scattering length corresponds to a weak attractive interaction between the 
atoms. This case is known as the BCS (Bardeen-Cooper-SchriefFer) limit. 
On the other side of the resonance the scattering length is positive. In the 
EEC (Bosc-Einstcin condensation) limit the interaction is strongly attrac- 
tive and the fermions form deeply bound molecules. For this reason the 
unitarity limit a — > oo is also known at the BCS/BEC crossover. 

The unitarity limit is of interest to QCD practitioners for a for a number 
of reasons: 

• The unitarity limit provides an approximate description of dilute 
neutron matter. The neutron-neutron scattering length is a„„ = 
— 18 fm, and the effective range is r„„ = 2.8 fm. This means that 
there is a range of densities, relevant to the outer layers of neutron 
stars, for which the interparticle spacing is large compared to the 
effective range, but small compared to the scattering length. 

• The Fermi gas at unitarity is a high Tc superconductor. There is an 
attractive interaction in the spin singlet channel which leads to s- 
wave superconductivity below some critical temperature Tc. In the 
unitarity limit the only energy scale in the problem is the Fermi 
energy Ep, and we must have ksTc = aEp with some numerical 
constant a. Quantum Monte Carlo calculations (and experimental 
results) indicate that a ~ 0.15 [1,2], much larger than in ordinary 
(or even high Tc) electronic superconductors, but comparable to 
what might be achieved in color superconducting quark matter [3] . 

• The limit a ^ oo corresponds to a non-relativistic conformal field 
theory [4]. In the unitarity limit there is no scale in the problem 
(other than the thermodynamic variables temperature and den- 
sity). Indeed, one can show that the theory is not only scale invari- 
ant, but invariant under the full conformal group. This raises the 
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question whether there are any physical consequences of conformal 
symmetry that go beyond results that follow from scale invariance. 
It also raises the possibility that a holographic description, similar 
to the AdS/CFT correspondence, can be obtained [5,6]. 
• Non-relativistic fermions at unitarity behave as a very good fluid 
and show interesting transport properties, including a very small 
shear viscosity. Kinetic theory suggests that the shear viscosity is 
inversely proportional to the scattering cross section, and reaches 
a minimum at unitarity. This expectation is confirmed by experi- 
ments that demonstrate large elliptic flow and a very small damping 
rate for collective oscillations [7,8]. 



2. Equilibrium Properties 

We begin by analyzing equilibrium properties of the dilute Fermi gas at 
unitarity. If the temperature is large, ksT > Ep, then the scattering cross 
section is regularized by the thermal wave length, and the effective inter- 
action is weak. Here the Fermi energy is defined by Ep = {3Tr'^nf^^/i2m), 
where n is the density, and m is the mass of the atoms. In the high temper- 
ature regime the equation of state is well described by the Virial expansion, 
and the system has single particle excitations with the quantum numbers 
of the fundamental fermions. In the regime A;bT ~ Ep the interactions 
are strong. As noted above, superfluidity occurs at ksTc 0.15Ep. Below 
the critical temperature the excitations are Goldstone bosons. In following 
section we will discuss the efi'ective theory of the Goldstone bosons, and 
relate the parameters in the effective lagrangian to static properties of the 
system. 



2.1. Low Energy Effective Theory and Density Functional 

The Goldstone boson field can be defined as the phase of the difermion 
condensate {iptp) = e^^'^ \ {'iptp) \ ■ The effective Lagrangian at next-to-leading 
order (NLO) in derivatives of (fi and the external potential is [9] 

C = com^/^X^/^ + cim^^^^^ + ^ f(VV)' - QmVVl Vx , (2) 

X yfn L J 

where we have defined 

Here, /i is the chemical potential and t) is an external potential. The 
functional form of the effective lagrangian is fixed by the symmetries of 
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the problem, GaHlcan invariance, U{1) symmetry, and conformal symme- 
try. The NLO effective lagrangian is characterized by three dimensionless 
parameters, co,ci,C2. These parameters can be related to physical proper- 
ties of the system. The first parameter, cq, can be related to the equation 
of state. We have 

25/2 



^o = Y^^, (4) 

where ^ determines the chemical potential in units of the Fermi energy, 

= ^Ep. The two NLO parameters ci,C2 arc related to the momentum 
dependence of correlation functions. The phonon dispersion relation, for 
example, is given by 



2C [c^ 



tC2 



+ 0{qHog{q^)) 



where Vs 



'^/3vp and vp = kp/m. The static susceptibility 

involves a different linear combination of Ci and C2, 

mkp 



1 + 27r2y2C ( ci - lc2^ ^ + 0{q' log(g2)) 



(5) 



(6) 



(7) 



Higher derivative terms in the effective lagrangian can also be used to com- 
pute the energy of inhomogeneous matter. At NLO in an expansion in 
derivatives of the density we find the following energy density functional [10] 



£{x) = n{x)V{x) + 



3 ■ 2- 



55/3, 



mc. 



2/3 



n{x) 



5/3 _ 4 2ci + 9c2 {Wn{x)y 



45 mco 



n{x) 



(8) 



12 C2 



W^n{x) . 



5 mco 

The first two terms correspond to the local density approximation (LDA) 
and the terms proportional to ci and C2 are the leading correction to the 
LDA involving derivatives of the density. 



2.2. Epsilon Expansion 

At unitarity the determination of Ci and C2 is a non-perturbative problem, 
and we will perform the calculation using an expansion around d = A — e 
spatial dimensions [11,12]. Our starting point is the lagrangian 



V2 

«9o + 0-3 — 
2m 



+ /U^Vs* + Ni^a+^(f) + h.c.) - 

Co 



(9) 
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Fig. 1. Leading order contributions to tiie effective potential in the epsilon expansion. 
Solid lines denote fermions propagators, dashed lines denote boson propagators, and the 
cross is an insertion of the chemical potential. 



where ^ = (■(/'|, V'j)"^ is a two-component Nambu-Gorkov field, ai are Pauli 
matrices acting in the Nambu-Gorkov space, a± = (ai ± zcr2)/2, is a 
complex boson field, and Co is a coupling constant. In dimensional regular- 
ization the fermion-fermion scattering length becomes infinite for 1/Co ^ 0. 

The epsilon expansion is based on the observation that the fermion- 
fermion scattering amplitude near d = 4 dimensions is saturated by the 
propagator of a boson with mass 2m. The coupling of the boson to pairs of 
fermions is given by 



m 



27T 



e/4 



(10) 



In the supcrfluid phase 
the boson field as 6 = d 



acquires an expectation value cpo ~ {4>) . We write 
I -|- gip. The lagrangian is split into a free part 



«9o +0-3- h (?io(o-+ + (T_) 

zm 



Am 



and an interacting part Ci + Cct, where 

Ci = g (^V+^v- + h.c) + iJ.^^a3^ + , 

V2\ 



(11) 

(12) 
(13) 



Note that the leading self energy corrections to the boson propagator gen- 
erated by the interaction term Cj cancel against the counterterms in Cct- 
The chemical potential term for the fermions is included in £j rather than 
in jCq. This is motivated by the fact that near d = 4 the system reduces to a 
non-interacting Bose gas and /i ^ 0. We will coimt fi as a quantity of 0{e). 
The Feynman rules are quite simple. The fermion and boson propagators 
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Fig. 2. Scalar self energy at LO in the epsilon expansion. 
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oil ' 



(14) 
(15) 



where E"^ = Cp + fPo ^^'i^ = I (2m). The fermion-boson vertices are iga"^ 
Insertions of the chemical potential are i/icrs. Both g"^ and /x are corrections 
of order e. 

In order to determine co , ci , C2 we have to compute three physical observ- 

ablcs. Wc have studied ^ = ji/Ep, and the curvature terms in the phonon 
dispersion relation and the static susceptibility. The universal parameter ^ 
was originally calculated by Nishida and Son. They computed the effective 
potential to NLO in the epsilou expansion, see Fig. 1. The derivative of 
the effective potential with respect to /i determines the density n, and the 
relation between n and /i fixes ^. The result is 



3/2 



l + ielog(e)-i(12C- 



5 + 51og(2))e + 0(e2) 



(16) 



with C = 0.144. The phonon dispersion relation can be extracted from the 
scalar propagator. Wc introduce a two-component scalar field $ = (ip, (^*). 
The scalar propagator satisfies a Dyson-Schwinger equation [13] 
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n 



(17) 



At LO in the epsilon expansion the self energy is determined by the di- 
agrams shown in Fig. 2. NLO contributions were calculated in [10]. The 
phonon dispersion relation is 



Po 



1 + ^ 



1 



1-1 
4 



(18) 
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Fig. 3. Leading order contributions to the static susceptibility. The wiggly hne denotes 
an external current. The double line is the scalar propagator defined in equ. (17). 

We note that the dispersion relation curves up (unhke ^He, but similar to 
weakly interacting Bose gases). This implies that there is — > + decay. 
Finally, we can determine the static susceptibility. Computing the diagrams 
in Fig. 3 we get [10,14] 




The coefficient co follows from the result for ^ = 0.475 at NLO in the 
e-expansion) using equ. (4). Matching equ. (18,19) against equ. (5,7) gives 
C2 = and ci/cq = 3/8 — e/4. The corresponding energy density functional 
was studied in [10]. Compared to a free Fermi gas the local density term is 
reduced by a factor ~ 2 (the interaction is attractive), while the gradient 
correction proportional to (Vn)^/n is enhanced by a factor ^ 2. 

3. Transport Properties 

In the following we will discuss transport properties of the Fermi gas at 
unitarity. The interest in non-equilibrium properties arises from the obser- 
vation that transport coefficients are much more sensitive to the strength 
of the interaction than thermodynamic quantities. A renewed interest in 
transport properties was also sparked the AdS/CFT correspondence and 
the experimental limits on the shear viscosity of the quark gluon plasma 
obtained at RHIC. In the following we shall focus on the shear viscosity of 
the Fermi gas at unitarity. Close to equilibrium the (coarse grained) energy 



August 5, 2008 22:23 



WSPC - Proceedings Trim Size: 9in x 6in caQCDOS 



8 




Fig. 4. Leading order processes that contribute to the shear viscosity at low temperature 
(Fig. a) and high temperature (Fig. b). Dashed lines are phonon propagators and solid 
lines are fermion propagators. 

momentum tensor can be written as 

Tij = {P + e)viVj - PSij + 6Tij , (20) 

6Tij = -viViVj + ViVj - ^SijV •-?) + ••-, 

where e and P are the energy density and pressure, and Vi is the local flow 
velocity. The first term is the ideal gas contribution, and STij is the leading 

order (in gradients of Vi) dissipative correction. The traceless part of ST^j 
is proportional to the shear viscosity rj. 



Tii=< I (21) 



3.1. Kinetic Theory 

We first consider the case that the fluid is composed of weakly interacting 
quasi-particles. In the unitarity limited Fermi gas this is the case at T <C Tc 
(phonons) and T ^ Tc (atoms). In these limits we can compute the shear 
viscosity using kinetic theory. In the following we will concentrate on the 
low temperature case discussed in [15]. In kinetic theory the stress-energy 
tensor is given by 

d^P PiPj 
(2^)3 Ep 

where fp is the distribution function of the phonons, Vg is the speed of sound. 
Pi is the momentum and Ep the quasi-particle energy. Close to equilibrium 

fp — fp''' + ^fpj where fp^^ is the Bose- Einstein distribution and Sfp is a 

small departure from equilibrium. Wc write Sfp = —x{p)fp'\^ + fp'^)/T- 

In the case of shear viscosity wc can further decompose 

1 2 ^ 

X{P) = 9{p){PiPj - 3^yP^)(ViUj + VjVi - -6ijV ■ v) . (22) 

Inserting equ. (22) into equ. (21) we get 
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The non-equilibrium distribution g{p) is determined by the Boltzmann 
equation 

= f „-.V/, = C[U (24) 

relating the rate of change of the distribution function fp to the collision 
operator C[fp]. The 2 <-> 2 collision integral is given by 

r \f] = J- [ ^'fe ^'fe' '^'p' (2^) 

^^^^■'^^ 2eJ {2TTf2Ek {2^f2Ek> {2T:f2Ep. ^ ' 

X {2TTf 6^^\p +k~p' - k')\M\''D2^2., 

where Z?2v^2 contains the distribution functions and \Jvl\ is the 2^2 
scattering amplitude shown in Fig. 4. The three and four-phonon vertices 
are fixed by the effective lagrangian (2). Linearizing Z)2«2 in 5fp one finds 

i)2«2 = ^/l?Vr(l + /r)(l + /r) i.x{p)+x{k)-x{p')-x{k')) . (26) 

There are a variety of methods for solving the linearized Boltzmann equa- 
tion. A standard technique is based on expanding g{jp) in a complete set 
of functions. A nice feature of this method is that the truncated expansion 
gives a variational estimate 

where hg is a set of expansion coefficients, Aqq is a normalization integral, 
and Mst arc matrix elements of the linearized collision operator. For the 
best trial function we find [15] 

,7/5 = 7.7x10-*'^^^, (28) 

where ^ is the universal parameter introduced in Sect. 2.1 and we have nor- 
malized the result to the entropy density s of a weakly interacting phonon 
gas. A similar estimate can be obtained in the high temperature limit. 
In this case the relevant degrees of freedom are atoms, and the dominant 
scattering process is shown in Fig. 4b. The result is [16,17] 

The high and low temperature limits of ry/s are shown in Fig. 5, together 
with the proposed lower bound 77/s = l/(47r) [18] and experimental data 
which we will discuss in the next section. 
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Fig. 5. Viscosity to entropy density ratio of a cold atomic gas in the unitarity limit. 
This plot is based on the damping data published in [8] and the thermodynamic data 
in [19,20]. The dashed line shows the conjectured viscosity bound »)/s = l/{An), and the 
solid lines show the high and low temperature limits. 



3.2. Hydrodynamics 

Hydrodynamics describes the evolution of long-wavelength, slow- 
frequency modes. The hydrodynamic description remains valid even if there 
is no underlying kinetic theory. The hydrodynamic equations follow from 
conservation of mass (particle number), energy and momentum. In a non- 
relativistic system the equations of continuity and of momentum conserva- 
tion are given by 

(30) 

mn— + mn yv ■ V j v = —VP — riS/V, (31) 

where n is the number density, m is the mass of the atoms, v is the fluid 
velocity, P is the pressure and V is the external potential. In an ideal fluid 
the equation of energy conservation can be rewritten as conservation of 
entropy, 

il-+V.(n,st/) = 0. (32) 
at 

A non-zero shear viscosity leads to dissipation, converting kinetic energy to 
heat and increasing the entropy. The shear viscosity of the dilute Fermi gas 
in the unitarity limit can be measured by studying the damping of collective 
modes in trapped systems [21]. The frequency of these modes agrees well 
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the prediction of ideal hydrodynamics. The dissipated energy is given by 



2 



E = J (fxf]{x) ^diVj + djVi - ^5^j^kVkj ■ (33) 

The damping rate is given by the ratio of the energy dissipated to the total 
energy of the collective mode. The kinetic energy is 

Ekin = Y / d^xn{x)v\ (34) 

If the damping rate is small both E and Ekin can be computed using the 

solution of ideal hydrodynamics. We recently performed an analysis [22] 
which is based on measurements of the damping rate of the lowest radial 
breathing mode performed by the Duke group [8]. We showed that can 
relate the dimensionless ratio T/ui, where F is the damping rate and u) is 
the trap frequency, to the shear viscosity to entropy density ratio 

..i,.<3.,. (-)(£-)(-). 

Here N is the total number of particles in the trap (2 • 10^ in [8]), ^ is the 
universal parameter defined in Sec. 2.1, E/Et=o is the ratio of the total en- 
ergy to the energy at T = (which can be extracted using a Virial theorem 
from the measured cloud size), and S/N is the entropy per particle (which 
is measured using adiabatic sweeps to the BCS limit [20]). The results are 
compared to theoretical prediction in the high and low temperature limit 
in Fig. 5. The data show a minimum near T/Tp ~ 0.2. At the minimum 
ri/s ^ 1/2. This should probably be considered as an upper bound, since 
dissipative mechanism other than shear viscosity may be present. In the 
high T limit there is fairly good agreement with kinetic theory. The tem- 
perature dependence implied by the low T prediction is not seen in the 
data. This is maybe not very surprising, since the mean free path in the 
low T regime quickly exceeds the system size. 



4. Outlook 

There are many promising directions for further study. Clearly, it is desir- 
able to obtain additional experimental constraints on the shear viscosity, 
and to improve the theoretical analysis of the existing data sets. It would 
also be interesting to confirm that the bulk viscosity vanishes in the nor- 
mal phase, and to measure the thermal conductivity. We would also like 
to improve the theoretical tools for computing transport properties in the 
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interesting regime near Tc- There are some recent ideas for applying holog- 
raphy and the AdS/CFT correspondence to Galilean invariant conformal 
field theories [5,6], but there are also many purely field theoretic methods 
(e expansions, large N methods) that have yet to be pursued. 
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